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ABSTRACT 

The majority of atomic hydrogen Br7 emission detected in the spectra of 
young stellar objects (YSOs) is believed to arise from the recombination regions 
associated with the magnetospheric accretion of circumstellar disk material onto 
the forming star. In this paper, we present the results of a K-band IFU spectro- 
scopic study of Br7 emission in eight young protostars: CW Tau, DG Tau, Haro 



1 Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the Particle Physics and Astronomy 
Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), CNPq (Brazil), and CONICET (Argentina). 
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6-10, HL Tau, HV Tau C, RW Aur, T Tau and XZ Tau. We spatially resolve Br7 
emission structures in half of these young stars and find that most of the extended 
emission is consistent with the location and velocities of the known Herbig-Haro 
flows associated with these systems. At some velocities through the Br7 line pro- 
file, the spatially extended emission comprises 20% or more of the integrated flux 
in that spectral channel. However, the total spatially extended Br7 is typically 
less than ~10% of the flux integrated over the full emission profile. For DG Tau 
and Haro 6-10 S, we estimate the mass outflow rate using simple assumptions 
about the hydrogen emission region, and compare this to the derived mass accre- 
tion rate. We detect extended Br7 in the vicinity of the more obscured targets in 
our sample and conclude that spatially extended Br7 emission may exist toward 
other stars, but unattenuated photospheric flux probably limits its detectability. 

Subject headings: stars: pre-main sequence — stars: winds, outflows — stars: 
formation — stars: individual (CW Tau, DG Tau, Haro 6-10, HL Tau, HV Tau 
C, RW Aur, T Tau, XZ Tau) 



Introduction 



H I emission lines are one of the defining characteristics of the classification of pre-main- 
sequence Sun-like sources known at T Tauri stars (TTS). Still in the midst of formation, the 
less evolved TTS, known as classical TTS (cTTS), are surrounded by optically thick disks of 
gas and dust. In most cases, these young pre-main-sequence stars are still interacting with 
and accreting matter from the innermost regions of their disks via stellar magnetic fields. 
In this magnetospheric accretion paradigm, the stellar magnetosphere guides disk material 
from the inner disk onto the stellar surface through magnetic channels. The gas travels 
along these channels or so-called accretion columns near free-fall velocities, terminating in 
an accretion shock at the stellar surface. It is generally accepted that the gas is heated and 
ionized prior to and after reaching the stellar surface and that the characteristic H I emission 
lines resul t, in part, from recombining and accreting hydrogen gas confined to these m agnetic 
channels jLvnden-Bell & Pringle!ll974i llJchida k Shibatalll984 : iBertout et allll988h . 



Balmer Ha emission is the dominant H I feature present in the optical spectra of cTTS, 
and the emission li ne strength (or line wid th) is often invoked as a measure of accretion rates 
for these sources (IMuzerolle et al.lll998al ). Spatially resolved observations of Ha emission 
lines in these young stellar objects (YSOs) also show it to be a strong component of the 
optical line emission from outflows, indicating that the underlying stimulation mechanisms 
for the H I lines is likely to be a combination of phenomena. Magnetospheric accretion models 
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can successfully reproduce many aspects of the H I emission features detected in the spectra 
of young st ars. However, they nota bly fail to account for the highest velocity gas in the H I 
line wings f lMuzerolle et al.l Il998a| ). Spectro-astrometric observations of the Pa/3 emission 
feature in the cTTS DG Tau show that the high velocity blue-shifted gas (v > -200 km s~ x ) 
is extended in the same direction as [Fe n] at 1.644 /xm, a forbidden emission line feature that 
traces the known outflow in the DG Tau system. This illustrates why a model producing 
H I emissio n from accretion fu nnels does not account for high velocity gas forming the H I 
line wings ( jWhelan et al.ll2004l ). While the high velocity gas is spatially shifted by as much 
as 0."5 around DG Tau (or ~ 70 AU at the distance of Taurus- Auriga) , the majority of the 
Pa/3 emission remains coincident with the source. Although these results demonstrate that 
there are multiple processes for stimulating H I emission lines, the bulk of the emitting gas 
does seem to arise from radii within 14 AU of the central source constraining the emission to 
the magnetospheric accretion columns, inner gaseous disk, and the base of disk winds and 
outflows. 

In the infrared, Br7 (2.16 /xm) e mission serves as a surrogate for Ha as a signpost for 
circumstellar disk accretion in TTS (INajita et al.lll996[ ). Br7 line luminosities appear to 
correlate with mass accretion luminosity in brown dwarfs, cTTS, and Herbig Ae/Be stars 
fHAEBEs: [Muzerolle et aDll998al [Natta et al.ll2004[ llVlohantv et al.ll2003l bood ). Moreover, 
Br7 and other lines in the infrared are less affected by optical depth effects that have proven 
problematic for using Balmer series lines to infer temperatures, densities, and geometries of 
the emitting gas. Historically, the problem of predicting the spectra emergi ng: from a recom- 
bining hydrogen gas has been divided in to two distinct cases, A and B ( IBaker fe Menzel 
19381 ). Case A theory applies to very low density gases that are optically thin to all tran- 
sitions of the hydrogen atom, including the ultraviolet photons associated with the Lyman 
series transitions. Case B theory, which applies to higher density gases which are optically 
thick to UV photons but optically thin to all n > 2 transitions, is often applied to the envi- 
ronments of T Tauri stars. In a Case B model of a recombining atomic hydrogen gas, Br7 
is ~0.8-l% of the flux of Ha for a wide range of densities and temperatures (10 2 < n e < 
10 6 cm -3 , 5000K < T < 20 000K). Ha is know n to be a strong component in optical line 



emission from YSO outflows ( INajita et al.lll996f ). it seems natural that a corresponding com- 
ponent of the Br7 emission would also arise from the outflows. Yet, to date, there has been 
little evidence in the literature for spatially resolved Br7 emission in the vicinity of YSOs, 
and hence nearly all Br7 emission is assumed to arise from magnetospheric processes within 
the inner accretion zone. In fact, infrared interferometric observations of HAeBe stars reveal 
that the Br7 often aris es from very compact locations with in the dust sublimation radius of 
the circumstellar disk (jEisner et al.ll2009l ; iKraus et al.ll2008l ). Though, a small extended Br7 
emission component can not be ruled out based on these observations. Further interferomet- 
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ric and spectro-astrometric programs that seek to reveal the inner ~1 AU environments show 



that the compact B r7 is not always well modeled by disk emission alone (jEisner et al.ll2010 



Malbet et al.ll2010l ). The analysis implies that a non-negligible component from outflowing 
gas needs to be incorporated into the models to explain the Br7 emission structure. 

As an extension of the spectro-astrometric and interferometric studies mentioned above, 
imaging spectroscopy of these TTS can provide us with considerable insight into the geo- 
metric distribution of H I emitting gas in accreting systems. Three-Dimensional imaging 
spectroscopy techniques can help to disentangle the respective contributions to the H I emis- 
sion features. With just one pointing of a telescope, imaging spectroscopy with integral 
field units (IFUs) can provide three-dimensional x, y, A datacubes at high spatial resolution 
with simultaneous coverage of many emission lines of interest. There has recently been an 
increase in the capabilities for a daptive optics (AO) fed near IR integral field spectroscopy 



at 8- 10 meter class observatories (jEisenhauer et al.ll2000t iMcGregor et al.ll2003l ; lLarkin et al. 



20061 ). IFUs optimized for AO spectroscopy have the power to spatially resolve emission line 
structures with less than 0."1 extents over the full wavelength ranges sampled by typical 
IR spectrographs. As such, the new generation of IFUs provides the means to study the 
accretion and outflow environments in cTTSs. 

In this paper, we present detections of spatially resolved Br7 emission in YSO envi- 
ronments from data acquired using the Near IR Integral Field Spectrograph at the Gemini 
North Observatory. We report on Br7 arising from eight classical TTS systems, and par- 
ticularly highlight the spatially extended emission detected in four of these: DG Tau, Haro 
6-10 (also known as GV Tau), HL Tau and HV Tau C. 



2. Observations 

Observations of the eight CTTSs listed in Table 1 were obtained using the Near IR 
Integral Field Spectrograph (NIFS) at the Frederick C. Gillette Gemini North Telescope 
on Mauna Kea, Hawaii. NIFS is an image slicing IFU fed by Gemini's Near IR adaptive 
optics system, Altair, that is used to obtain integral field spectroscopy at spatial resolutions 
of <0."1 with a spectral resolvi ng power of R~5300 at 2.2 /im (as measured from arc and 



sky lines; IMcGregor et al.l 120031 ) . The NIFS field has a spatial extent of 3" x 3", and the 
individual IFU pixels are 0. "1x0. "04 on the sky. Data were obtained at the standard K-band 
wavelength setting for a spectral range of 2.003-2.447 /im. All observations were acquired in 
natural seeing of better than 0."7 for excellent AO correction. 



The data sets for this study were acquired for commissioning and system verification 
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of NIFS in October 2005 and February 2006, GTO time in December 2006, and in queue 
mode in February 2007 (see Table 1). For each observation, a standard set of calibrations 
were acquired using the Gemini facility calibration unit, GCAL. The raw IFU frames were 
reduced into datacubes using the NIFS tasks in the Gemini IRAF packag^j]. 



Beck et all (120081 ) discuss these NIFS data on DG Tau, HL Tau, HV Tau C, RW Aur, 
T Tau and XZ Tau in the context of spatially resolved molecular hydrogen emission lines. 
Hence the observational details, calibration, and data reduction information is described in 
great detail in that paper and excluded here. To study the Br7 emission in the young stars, 
the absorption features from the AO stellar type telluric calibration stars were removed by 
fitting and dividing Voigt absorption profiles in the 2.16 /im spectral region and cleaning the 
calibration spectra for any small residuals. Observations of Haro 6-10 were acquired with 
the laser-fed AO system using the R~16.5 mag southern component in this l."2 binary as 
the laser tip-tilt reference star. These data were obtained in excellent laser-quality weather, 
photometric with better than 0."5 seeing. The CW Tau and H aro 6-10 data were processed in 



a similar manner to all other data, as described above and in lBeck et al.l (120081 ) . All or part 
of the data for each source in this project was observed during photometric conditions and 
were flux calibrated using K-band magnitudes estimated by comparison to the brightness 
of the AO standard star used for telluric correction. With the exception of Haro 6-10, the 
derived fluxes of the systems (combined in the case of multiples) were within 10-15% of 
published or 2MASS magnitudes. For Haro 6-10, the NIFS IFU spectra were compared to 
infrared K-band images acquired nearby in time for a complimentary project, the difference 
in flux was less than 0.2 magnitudes between the two observations. We estimate that our 
overall data flux calibration is good to ±10-15%. The data cubes were interpolated onto a 
square pixel grid with 0."05 spatial sampling, and the velocity channel steps through the IFU 
cubes are ~29 km s~ x pixels in extent at 2.20 /xm (~56 km s" 1 2 pixel resolution). The final 
reduced, combined, telluric corrected, and flux-calibrated datacubes of the Br7 line emission 
for each target are discussed in detail in the following sections. 



3. Spatially Extended Br7 Emission from YSO Environments 



The data acquired for this project, described in detail in the preceding section, were 
obtained with the goal of studying th e K-band fe a tures of molecular hydrogen emission. The 
IFU spectral data were discussed by iBeck et al.l ( 120081 ) in this context. We never expected 
to detect spatially extended Br7 emission in any CTTS. This was a serendipitous discovery, 



Information on the Gemini IRAF package is available at 



http://www.gemini.edu/sciops/data/dataIRAFIndex.html 
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revealed as we stepped through the raw velocity cube at Br7 wavelengths in the Haro 6-10 S 
data. The clear detection of spatially extended emission from the Haro 6-10 S jet prompted 
us to take a closer look at the Br7 emission in all TTS for which K-band IFU spectra 
had been obtained. As presented and discussed in the following, we have found significant 
spatially extended Br7 emission in four of the eight stars presented here: DG Tau, Haro 6- 
10, HL Tau and HV Tau C. We do not find appreciable spatially extended Br7 in CW Tau, 
T Tau, XZ Tau, or RW Aur. 

Figures 1 through 4 show the images of the spatially extended Br7 emission from DG Tau 
(Figure 1), Haro 6-10 (Figure 2), HL Tau (Figure 3) and HV Tau C (Figure 4). The 
panels in these figures show: a) - the continuum emission with contours of [Fe II] emission 
overplotted to demonstrate the outflow position and geometry, b) - continuum subtracted Br7 
line emission maps with contours of the continuum overplotted, c) - "point-source subtracted" 
spatially extended maps of Br7 emission with the contours of the continuum-subtracted Br7 
overplotted, and d) - images of the integrated blue-shifted Br7 emission only, with contours 
of the point-source subtracted Br7 emission overplotted. At the right side of each image is 
the key correlating the image intensity to flux. The continuum images in panels (a) were 
constructed by fitting a straight line to the continuum around the Br7 emission feature 
and integrating the linear fit through the velocity channels that correspond to the Br7 
emission. The images of the total Br7 emission presented in the b) panels were constructed 
by subtracting the linear fit to the continuum from the datacubes, then integrating in velocity 
over the Br7 emission feature. The "point-source subtracted" Br7 images shown in the c) 
panels were derived by normalizing the continuum (PSF) image to the peak flux in each 
velocity channel through the Br7 emission datacube, subtracting this scaled continuum image 
off of the Br7 cube, and integrating over the velocity extent to form an image of only the 
extended emission. The image of spatially extended blue-shifted Br7 emission shown in panel 
d) was made by integrating the 'point-source subtracted' image in 2-3 velocity channels of 
blue-shifted emission only. The d) panels show that the spatially extended Br 7 emission 
is stronger in the blue-shifted velocity channels. Because of Poisson statistics associated 
with the subtraction process, detection of spatially extended Br7 emission is less robust 
at distances of <0."1 from the central point-source. The Br7 emission in the vicinity of 
HV Tau C is quite weak, and all of the detected line emission is spatially extended. As a 
result, Figure 4 presents only the (a) and (b) panels for HV Tau C. 

For DG Tau and Haro 6-10 S, the majority of the Br7 emission that we detect, integrated 
over wavelength, is consistent with the stellar point-source image. For HL Tau and HV Tau 
C, the majority of the Br7 is not coincident with the central TTS. The integrated Br7 
emission from DG Tau follows the point-source continuum contours with little deviation 
(Figure lb). The spatially extended Br7 emission from DG Tau is detected in panel 1(c), 
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and seen clearly in the blue-shifted emission shown in panel 1(d). The Blue-shifted emission 
(Figure Id) extends to the south-west of DG Tau, at an orientation and velocity consistent 
with the known, collimated blue-shifted jet. Haro 6-10 also shows the bulk of the Br7 
emission arising from the two stellar point-sources. However, strong emission comes from 
the same location as the [Fe n] outflow emission from Haro 6-10 S (Figure 2c). In fact, for 
Haro 6-10 S, the Br7 arising from the outflow is ~20% of the total spatially integrated Br7 
emission in some blue-shifted velocity channels (Figure 2d). No spatially extended Br7 is 
detected toward the Haro 6-10 N component. 

Data on HL Tau were acquired with the 0."2 occulting disk blocking the central stellar 
point-source, however an appreciable amount of the detected Br7 emission is seen to deviate 
from the central stellar position. HL Tau shows strong spatially extended Br7 t hat follows the 



scatt e red light nebulosi ty revealed in the continuum emission (Figure 3b, 3c; iTakami et al. 



20071 ; Oose et al.lll997l ). The integrated blue-shifted Br7 from HL Tau (Figure 2d) shows 
very weak emission detected at a 4<r level of significance that corresponds precisely with the 
spatial location of the [Fe il l emission. HV Tau C is a system with a known circumstellar 
disk viewed nearly edge-on (IStapelfeldt et al.ll2003l ) and the stellar continuum flux is much 
fainter than the other sources. Curiously, we find no significant Br7 emission associated with 
the locations of the continuum flux. The Br7 emission from HV Tau C is quite weak. It 
is spatially extended from the continuum and it appears to only follow the location of the 
outflow seen in [Fe il] (Figure 4). 

The spatially extended Br7 emission revealed in Figures 1, 2, and 4 for DG Tau, Haro 6- 
10 S, and HV Tau C, lies precisely along the outflow axis of the known Herbig-Haro energy 
flows associated with these young stars. For HL Tau, the detected (low signal-to-noise) Br7 
emission in the blue-shifted component of the e mission (Figure 3d) arises from the same 
spatial location as the known blue-shifted outflow ( ITakami et al.ll2007l ). However in HL Tau 
the majority of the spatially extended Br7 emission (Figure 3c) is not appreciably shifted 
in velocity from the nominal stellar radial velocity; it is detected at much higher S/N, and 
appears to arise from Br7 emissi on from the central point-source that h as been scattered off 
of the wall of the outflow cavity fjClose et al.l 119971 ; ITakami et al.ll2007l ). 



Figure 5 shows the velocity profiles of the Br7 flux associated with the central point 
source for DG Tau, Haro 6-10, HL Tau and HV Tau C at the location of the peak con- 
tinuum emission (upper panel) and the profiles of Br7 emission extracted in 0."2 diameter 
apertures at "Position B" and "Position C" as designated in Figures 1 through 4 for each 
star. The spatially extended Br7 from DG Tau and Haro 6-10 S is blue-shifted in velocity by 
> 100 km s _1 with respect to the central point-source flux. HL Tau shows strong emission 
at the velocity of the point-source flux in "Position B," and very slight ~3-4<r detection of 
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flux from blue-shifted (~-200 km s 1 ) emission from the outflow. 

In DG Tau, Haro 6-10 S and HL Tau, the spatially extended Br7 emission that we detect 
from the outflows corresponds to the blue-shifted regions of the jets. This is consistent with 
the fact that blue-shifted jet components are flowing into our line of sight, and are hence less 
obscured by intervening circumstellar disk material. HV Tau C is the only source where both 
blue and red-shifted Br7 is detected from opposite sides of the outflow with respect to the 
(estimated) position of the central star. The extended Br7 from HV Tau C is quite weak in 
some regions, it is detected at a low signal-to-noise but follows the location of the extended 
outflow. Curiously, we find that the blue-shifted outflow lies t o the northeast, and th e red- 



shifted emission is to the southwest, which (as predicted by IStapelfeldt et al.l 120031 ) is at 
odds with the relative brightnesses of the lobes of the scattered light edge-on disk reflection 
nebula associated with this source. The southwestern lobe of the HV Tau C scattered light 
nebula is brighter, and was thus thought to be associated with the blue-shifted (closer) lobe 
of the outflow. Our data reveal that this is not the case, and the north-eastern, fainter lobe 
seems to be associated with more blue-shifted outflow emission. 

Figure 6 presents the standard deviation of continuum subtracted Br7 flux with in- 
creasing distance from the central star (plotted as a solid line) for HL Tau (a), DG Tau 
(b), Haro 6-10 (c) and HV Tau C (d). These curves were derived for each source by cal- 
culating the mean flux within 0."08 wide pixel annuli from the central stellar position and 
computing the standard deviation from the mean. This standard deviation of continuum 
subtracted line emission provides a measure of the uncertainty in the subtraction process, 
and thus can be used to gauge the S/N of the spatially extended line emission. Though, the 
continuum subtracted plots do have the extended Br7 flux included within the field. As a 
result, the measured standard deviations in the spatial bins with appreciable Br7 flux are 
correspondingly higher, and the accuracy of the subtraction process is thus under-estimated. 
Overplotted in each of the four panels is a dashed line that shows the magnitude of the 
point-source subtracted Br7 flux with increasing distance from the central stellar position. 
Comparison of the dashed curve with the solid curve provides an estimate of the measured 
S/N of the extended Br7 emission. The peak S/N for extended Br7 is ~7 for HL Tau, ~5 
for DG Tau, -25 for Haro 6-10 S and ~4 for HV Tau C. 

Figure 7 plots the same standard deviation of continuum subtracted Br7 flux with 
increasing distance from CW Tau, T Tau, XZ Tau, and RW Aur, and a dashed line is also 
over-plotted that presents the peak magnitude of the point-source subtracted Br7 emission 
from the central star (in the case of the multiple systems, the central star is assumed to be 
the brightest stellar component). The stellar companion position locations are apparent in 
the plots for XZ Tau and T Tau, while RW Aur's companion has a slightly greater separation 
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than presented in Figure 7d. For all four sources presented in Figure 7, no excess was seen 
in spatially extended Br7 emission at the location of the known Herbig-Haro outflows. For 
the cases of CW Tau, XZ Tau and RW Aur, no strong evidence of extended Br7 emission 
beyond a S/N of ~3 is found. Curiously, XZ Tau B showed weak Br7 emission associated 
with the position of the star, but XZ Tau A had no detectable Br7. This result is seemingly 
at odds with the proposition that XZ Tau A is the more actively accre ting star and the 



main driving source of the Herbig-Haro flow associated with this system (IKrist et al.l 120081 ) 
Although RW Aur exhibits very strong, centralized Br7 emission flux, RW Aur B showed no 
measurable Br7 emission associated with stellar mass accretion from its circumstellar disk. 

The central region of Br7 emission for T Tau N was saturated in the data, making a 
proper measurement of the continuum subtracted Br7 flux and point-source subtracted flux 
difficult. The surrounding spatial and spectral regions nearby are not saturated, and were 
used to estimate the 2.16 /im flux level based on the shape of the PSF. A cursory fitting 
analysis done at velocities on and off of the emission line using the A spectral type telluric 
calibrator as a PSF reference showed that the Br7 emission associated with T Tau N is not 
appreciably extended compared to the point-source continuum emission. The line emission 
is asymmetric around the T Tau S PSF, and is brighter to the northwest at the location of 
T Tau Sb. Hence, the Br7 line emission associated with the nearby 0."1 T Tau S a+b binary 
seems to arise preferentially from the b component, and it is stronger than the continuum flux 
(i.e., the T Tau Sb / Sa flux ratio is greater in Br7). This causes the apparent enhancement 
in Br7 emission at the position of the companion in Figure 7b. Extraction of the individual 
spectra of the blended components was not done because of the saturation of T Tau N, which 
would need to serve as a PSF calibrator. 



4. Br7 Estimates of YSO Mass Accretion and Mass Outflow Rates 



Muzerolle et al.l ( I1998bl ) showed that the Br7 line luminosity from young stars correlates 
with the mass accretion rate, as determined from UV excess emission. Hence, we now use 
our detected Br7 line fluxes to derive mass accretion rate for the observed targets. Table 2 
presents the total Br7 emission line fluxes (i.e., integrat ed over all velocity cha nnels) for 
each of the young stars in this study. The relation from iMuzerolle et al.l ( 1l998bl ) assumed 
all detected Br7 emission was associated with the stellar point-sources, so the line fluxes 
that we have included in Table 2 are the total integrated flux values, including emission 
detected in the outflowing gas. Also included in Table 2 are the adopted stellar parameters 
used to derive t he mass accretion rates: the stellar mass, temperature, luminosity and vi- 



sual extinction (IKenyon fc Hartmannlll995t White fc Ghezll2001t lHartigan fc Kenyonl 12003 
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Doppmann et al.ll2008l ). The stellar parameters seem consistent with our K-band spectra, 
so we do not rederive them from our data. 

XZ Tau B, T Tau S and Haro 6-10 N are "infrared luminous companions" (IRCs) 
to their respective primarie s , and the stellar pa rameters for these sources are much less 
certain (IKoresko et al.l 119971 ; IWhite &: Gheg|200ll ). The accretion activity and line of sight 



visual extinction toward the IRCs co uld also be variable (IGhez et al.lll99ll : iBeck et al.ll2001 



Leinert et al. 


2001; 


Beck et al. 


2004) 



the Br7 emi ssion that we detect arises not from the IRC but fr om the ~M-type T Tau Sb 



companion (IDuchene et al.ll2005l ). Haro 6-10 N was found by iDoppmann et al.l ( 120081 ) to 
have weak evidence of photospheric Na absorption features at 2.20 /im, with a late spectral 
type and an infrared veiling value estimated to be in the range of 12-15. Similarly, XZ Tau B 
also has strong infrared veiling and a poorly determined spectral type. As a result, we do 
not try to estimate mass accretion rates for T Tau S, Haro 6-10 N or XZ Tau B. Column 7 
of Table 2 prese nts the mass accretion r ates (M acc ) derived for all of the other stars, using 
the relation from iMuzerolle et al.l ( 11998b! ) for th e Br7 line luminosity to accretion luminosity 
and the Virial Theorem treatment put forth by lGullbring et al.l ( 119981 ). 



The mass accretion rates derived for the eight stars in this study that have well de- 
termined stellar parameters lie in the range from less than 4xl0 -10 M for RW Aur B, to 
1.5 xl0~ 7 M for T Tau N, with most sources in the range of 10 8 to 10 7 M . These 
values compare lie within the overall range of mass accretion rates derived for cTTSs 
( IMuzerolle et al.lll998bl ; iGullbring et al.lll998f ). For the most part, the mass accretion rates 
that we derive here are similar to M acc values for these stars that have been derived from 
previous studies, within the associated uncertainties (e.g., Hartigan et al. 1995; Muzerolle 
et al. 1998; White et al. 2004). However, direct comparisons between mass accretion values 
for DG Tau, particularly, show a large discrepancy in our study. DG Tau has in the past 
exhibited a mass accretion value on the high side of the range for TTSs; around or just under 
10 -6 M yr _1 . Our derived ma ss accretion rate is an order of magnitude le ss than many 



previously published estimates ( iHartigan et al.l Il995t IMuzerolle et al. 



998b). DG Tau is 



known to vary in flux and spectral characteri stics on short time scales (IBiscaya et al.lll997l ; 
Hessman Sz Guentherl Il997l ; iBary et al.l 120081 ) . Bary et al. (2008) found a nightly average 
Br7 equivalent width of 6.4 A for seven nights of observations with a standard deviation of 
20%. Therefore, the discrepancy we see in mass accretion rate toward DG Tau might be a 
result of both intrinsic stellar and accretion variability. 

Also presented in column 8 of Table 2 is the Br7 flux from the spatially extended outflows 
from DG Tau and Haro 6-10 S. For our analysis of Br7 emission seen in the outflows, we 
extracted a region of Br7 flux from the spatially extended jets from DG Tau and Haro 6- 
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10 (Figures Id and 2d). Multi-epoch observations have shown that the average transverse 
proper mot ion of an outflow is ~0."2 to 0."3 per year, or 28-42 AU at the ~140pc distance 
of Taurus ( jTorres et all 120091). Both the DG Tau and Haro 6-10 S outfl ows are viewed 
at 60-70° inclination angles jMovsessian fc Magakian1ll999l : IPvo et al.l l2003h . and the 28-42 
AU extent in the outflow is also roughly consistent with the corresponding annual motion 
of an average flow velocity of ~200 km s" 1 ; e.g., material in the flow would move 28- 
42AU/yr at the estimated geometry and velocity. For this reason, to capture the emitting 
flux from approximately one year of jet motion, we have extracted a segment, 0."3 in extent, 
of the spatially resolved Br7 emission along the jet axes for the DG Tau and Haro 6-10 
outflows. The width of the extraction box corresponds to a 0."2 radius encircling the jet. 
The Br7 emitting volume, Vhi, roughly corresponds to this area surrounding the jet, and is 
approximated by a cylinder of 28 AU (0."2) radius and 42 AU (0."3) height. The "one year" 
integrated Br7 outflow emission fluxes for DG Tau and Haro 6-10 S extracted from these 
apertures are 6.8xl0~ 16 and 8.5 xlO -16 erg s _1 cm -2 , respectively. 



Under the assumption that the spatially extended Br7 emission seen from DG Tau and 
Haro 6-10 S arises from thermally excited emission from optically-thin, post-shock regions of 



the outflows and behave as a Case B recombining gas (IBaker fc Menzell Il938l ; iBrocklehurst 



197ll ; iHummer fc Storeylll987l ; IStorey fc Hummerlll995l ). we use a simple analysis to measure 
the m ass outflow rate from these stars. Using the H I emission coefficients from lOsterbrock 
( 119891 ). the expression for the detected Br7 flux can be described as: 



F Br7 = 1.2 x 10- 28 (N c 2 Vhi/D 2 ) (1) 

where D is the distance to the emitting region in centimeters, N e is the electron density, 
and Vhi is the volume of the emitting region (MKS units). The electron temperature is 
assumed to be ~10 4 K, which is a reasonable estimate for the inner regions of YSO outflows. 
The mass of the emitting hydrogen in the spatially resolved regions of the outflow can be 
estimated as Mhi = m p N e VHi- Merging this equation for the mass with the above equation 
for the flux gives: 



M ffl = 1-5 x 10- 13 (F Br7 V ffl ) 1/2 D (2) 

We can thus solve for the mass of emitting hydrogen gas using the parameters of our measured 
Br7 flux, the selected emission volume, and the assumed 140 pc distance to the s tars, which is 



based on distances derived toward TTS within the Taurus star forming complex ( Torres et al. 
20091 ). For DG Tau, an H I mass of 1.2xl0 -8 M Q is derived, and for Haro 6-10 S this value 



is 1.4x10 8 M m . The flux extraction volumes were chosen to be the same, so the difference 
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in the atomic hydrogen outflow masses between the two stars is determined by the difference 
in their detected Br7 flux values. Having chosen the flux extraction volumes to correspond 
to the average annual proper motions of these jets, these estimates approximate the mass 
outflow rates for DG Tau and Haro 6-10 S in solar masses per year. 

Based on this analysis, we derive average electron densities on the order of a few 
xlO 4 cm -3 for the DG Tau and Haro 6-10 S outflows. These densities are consistent with 
the large values of 10 4 to 10 6 cm -3 for N e that are often found within the inner one hun- 
dred AU regions of young star outflows derived by inspecting [F e Hi and other forbidden 
emission species (IHartigan et al.l Il995l ; iBacciotti fe Eisloffell Il999t lHartigan fe Morsel 120071 ; 
Hartigan fe Hillenbrand! 120091 ). These electron densities are also very consistent with past 



values found in th e inner DG Tau high velocity blue-shifted jet ( IBacciotti fe Eisloffell 11999 
Coffev et al.ll2008h . 



The hydrogen outflow rates of 1.2 xlO^ 8 M Q yr^ 1 for DG Tau and 1.4xl0~ 8 M Q yr" 1 
for Haro 6-10 S represent lower limits for the true mass outflow levels derived in this manner. 
These mass outflow rates are underestimates of the true mass flow because only the fraction 
of the gas that has recently been heated by the shock radiates in the emission species that 
is studied. We also only detect the gas in the hi gh velocity component of the outflow , 
not in lower velocity flow surrounding the jet axis (IBacciotti et al.ll2000l ; iPvo et al.ll2003h . 
Moreover, the derived Br7 flux in the outflowing volume was not corrected for any line of 
sight extinction effects which might further serve to increase the derived mass outflow rate, 
particularly in the obscured Class I star, Haro 6-10 S. Thus, we predict that the true M out 
values are greater than 1.2xlO _8 M Q yr _1 for DG Tau and more than 1.4xlO _8 M yr _1 for 
Haro 6-10 S. 



5. Discussion 

This study shows that not all H I Br7 emission from classical TTS arises from magne- 
tospheric accretion processes within a few radii from the central star. We detect spatially 
resolved Br7 from DG Tau, Haro 6-10 S, HL Tau and HV Tau C, which represents 50% of 
the TTS systems in our sample. We do not spatially resolve Br7 line emission in the envi- 
ronments of XZ Tau, RW Aur, T Tau and CW Tau. Two stars within these latter systems, 
RW Aur B and XZ Tau A, exhibit no Br7 emission at all. 

In two of the blue-shifted velocity channels, the spatially extended Br7 from Haro 6-10 S 
makes up ~20% of the total spatially integrated line flux. Integrated over the whole velocity 
width, about 10±2% of the Br7 emission from Haro 6-10 S is spatially resolved at distances 
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of greater than 0."1 (14 AU) in the extended outflow (entirely from the blue-shifted velocity 
component). The spatially extended Br7 emission makes up ~2% of the total line flux for 
DG Tau. All of the detected Br7 flux seems to arise from the outflow for HV Tau C, but 
the star is not seen directly because the continuum flux is measured only from the scattered 
light nebulosity. The inner magnetospheric accretion region of HV Tau C may be shielded 
by the inner rim of the central circumstellar dust disk in our edge-on viewing orientation, 
hence in this case we might not see the strong central Br7 component in the scattered light 
nebulosity. HL Tau has considerable Br7 emission scattered off of its surrounding outflow 
cavity walls, and the fraction of spatially extended Br7 emission is estimated to be 18±7% 
of the line flux from the central position. Observations of HL Tau were acquired with an 
occulting disk in the beam, so the total point source flux is estimated from the PSF shape 
in the target acquisition setup images and the resulting value is significantly less certain. 

The majority of the integrated Br7 flux that we measure is spatially unresolved from the 
position of the central stellar sources in our data, with the noted exception of the edge-on disk 
system HV Tau C. NIFS data can spatially resolve the bright Br7 emission beyond ~0."1 
from the star, or about ~14 AU. Additional Br7 emission likely arises from the outflows in 
regions closer than our resolution limit, where it cannot be detected. Overall, the extended 
Br7 line emission that we detect beyond ~14 AU distances from the parent stars comprises 
anywhere from a few percent (DG Tau) to all of the detected line flux from these systems 
(HV Tau C). 

All of the sources that we have studied here are known to drive Herbig-Haro outflows. 
It is not clear why some stars exhibit extended Br7 emission, while other sources with strong 
and collimated outflows do not. However, from information presented in Table 2, we see that 
the stars that have appreciable spatially extended Br7 emission also have stronger estimated 
levels of optical visual extinction, A v , toward the stellar photosphere. The stars where we 
do not find extended Br7 emission all have lower estimated optical obscurations. In the 
systems where extended Br7 emission is seen, the stronger levels of stellar continuum flux 
attenuation from the high visual extinction may make the weak spatially extended emission 
easier to detect. The stars that do not exhibit extended Br 7 emission have brighter stellar 
continuum flux, which may be a result of less obscuration by natal material because of a 
slightly older evolutionary state, or a perhaps a more inclined viewing geometry that directly 
reveals more of the central photosphere of the star. Thus, spatially extended Br7 emission 
may exist toward the other stars, but the bright continuum flux might prevent us from 
detecting it. 

Several recent spectro-astrometric and interferometric investigations of inner YSO disks 
have sought to spatially resolve the atomic hydrogen associated with the magnetospheric 
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recombination regions w ithin 1 AU from the target stars flWhelan et al.1 12004c lEisner et al. 



20091 ; iKraus et al.ll2008l ). The spectro-astrometric study of IWhelan et al.1 ( 120041 ) revealed 



that an appreciable amount of Pa/3 emission from the blue and red-shifted velocity profiles 
arises from spatially extended distances from a TTS. Three of their four targets exhibited 
spatially resolved Pa/3 emission, two revealed evidence for spatially extended bi-polar red and 
blue-shifted components. This is an important finding, since the broadened emission line 
wings seen in the H I profiles from cTTS have never been well explained by magnetospheric 
accretion models. Evidence from this Pa/3 study and now our Br7 project suggest that the 
H I emission from TTS does have non-negligible components from the outflow on spatial 
scales that were unresolvable prior to the current generation of sensitive instrumentation on 
8-10 meter class telescopes. 

Thus far, spectro-interferometric studies of Br7 that can spatially resolve the inner 
magnetospheric accretion reg i on hav e largely concentrat e d on the brighter H erbig Ae/Be 
stars (HAEBEs; lEisner et al.1 J2009I ); IKraus et al.l (bo08h : iMalbet et all j2007h ). Interest- 
ingly, most interferometric data reveals that the Br7 emission surrounding the HAEBEs 
arises from a centralized location that is spatially more compact than the infrared contin - 
uum emission, with perhaps some contribution from an outflowing wind (IMalbet et al.ll2007l ). 
The compact H I likely arises predominantly from the accretion-driven processes within the 
star-disk boun dary; either from t he central stellar mass accretion engin e or from a stellar 
or disk wind ( iKraus et al.l l2008t IMalbet et al.l 120071 : lEisner et al.l |2010| ) . Further work is 
warranted to better reveal the central emission location of H I from TTSs, and the fraction 
of H I emission that is truly spatially extended from the parent star in comparison to the 
higher mass HAEBEs. 

Although we do detect spatially extended Br7 emission from 50% of our sample targets, 
the bulk of the line emission from most stars does arise within ~14 AU from the central 
unresolved point sources. Whil e the spatially exten ded H I seems to affect the wings of 
the velocity line profile shapes (IWhelan et al.ll2004l ). it is likely that the majority of the 
low velo city H I is e mitte d via magnetospheric processes in the central disk-star accretion 
region. iBary et al.l ( 120081 ) presented a large-scale multi-epoch near-infrared spectroscopic 



survey of 15 actively accreting TTS in Taurus-Auriga. They compared the H I line ratios of 
16 Paschen and Brackett series emission features and found little variations in these ratios 
from epoch to epoch and from source to source. For the first time, a statistically signifi- 
cant correlation was found between the observed H I line ratios in TTS systems and those 
predicted by Case B line reco mbination theory (IBaker fc Menzell Il938l ; iHummer fc Storey 
19871 ; IStorey &; Hummer! Il995l ) for a tightly constrained gas temperature and electron den- 
sity. While the range of electron densities Bary et al. find for the emitting H I gas agrees 
with that predicted by magnetospheric accretion models, the temperatures are substantially 
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lower than expected for accreting gag^|. These results suggest that the emission is from a 
recombining gas not in local thermodynamic equilibrium and is instead being stimulated by 
a non-thermal process. In addition, these results may also suggest a combination of emission 
sources for the emitting H I gas in the central regions and that the observed H I line fluxes 
are a superposition of emitting gas with different physical conditions. As such, one might 
not necessarily expect a linear correlation between detected Br7 flux from the central star 
and the liberated accre tion luminosity measured from the UV excess emission, as found by 
Muzerolle et al.l (jl998bl ). However, the correlation does exist and the mass accretion rates we 
derived from our measured Br7 line fluxes are consistent with those found in past studies of 
TTS. Taken in combination, these results further confirm the connection between accretion 
and outflow activity. 

Deriving accurate mass outflow values is notoriously difficult. This is in part because the 
emission line regions in outflows measure only material in the post-shock area of a jet and do 
not sample outflow material ahead of the shock. Shocks from outflows and their associated 
optical and infrared emission lines are obvious observational manifestations of mass outflow 
activity, but they are the result of fast material plowing into slower clumps of gas in the 
outflow path. Hence, material can (and does) exist in the outflow without these detectable 
emission features if there is no slower material in the flow path to cause the shock. Methods 
described in the literature for measuring mass outflow rates from classical TTSs rely largely 
on measuring the mass flux from forbidden emission line speci es and converting th i s to an 



outflow rate based on estimates of the gas ionizat i on fraction ( Hartigan et al 



20021 ; iDougados et all 120021 ; iBacciotti et al.l 120021 ; ICoffev et al 



2008 



1995 



Cabrit 



Agra-Amboage et al. 



20091 ) . The different observational techniques derive mass outflow to mass accretion ratios in 
the range of 1 to 10%. In some cases, multi-wavelength or emission line studies of the same 
system find M out /M acc values that differ by factors of ~10, even for the same system (e.g., 
DG Tau, see below). Measurement of mass outflow rates using atomic hydrogen emission 
species has been difficult in the past because of the strong optical depth effects in optical 
Balmer emission species. The infrared Brackett features are less optically thick, and may 
serve as a more robust tracer of the atomic hydrogen gas mass in the young star outflows. 

In the previous section, we derived the first masses of ionized hydrogen using Br7 
emission for the DG Tau and Haro 6-10 S outflows. Comparison of the mass accretion 
rates with our derived atomic hydrogen mass outflow rates for DG Tau and Haro 6-10 S 
(see Table 2) reveals that the outflowing H I mass is estimated at ~10-15% of the mass 
accreting onto the central star. Moreover, our derived mass outflow rates are also lower limits, 



Barv et al.1 (|2008l) found the most likely temperature and density ranges of T < 2000 K and 



10 9 cm" 3 < n e < 10 1 
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implying that the actual ratio of M out /M acc may be larger than 10-15% for these two sources. 
These ratios are consistent with M out /M acc values predicted by the theory of mass accretion 
efficiency and magnetocentrifugal launching mechanisms (i.e., disk winds and accretion- 
drive n stellar winds) for young star outflows (|Koenig]||l99ll ; iMatt fc Pudritzll2005l ; ICranmer 
20081 ). and are on the high side in comparison with ratios found from other observational 



investigations (Edwards et al 



1987 



Hartigan et al. 



1995 



Cabrit 



2002 



Dougados et al. 2002 



Bacciotti &: Eisloffel 



1999 



Coffey et al. 



2008 



Agra-Amboage et al.ll2009i ). 



It is interesting that we find M out /M acc values for DG Tau and Haro 6-10 S that seem in 
the high range compared to ratios derived from other observations. Additionally, the M out 
rates we find are lower limits, as discussed in §5. However, the uncertainties in the derivation 
of both the M out and M acc values can be very large. Our lower limit for the value of M out 
in DG Tau and Haro 6-10 S is derived from straightforward assumptions on the physics of 
LTE atomic hydrogen in recombination regions. As noted above, the hydrogen recombination 
emission area very close to th e central star is po ssibly excited by non-thermal processes in the 
inner magnetospheric region (IBary et al.ll2008l ). Hence, it seems feasible that these non-LTE 
conditions may extend into the inner regions of the outflow, making our H il-like analysis of 
hydrogen emitting mass correspondingly uncertain. This is likely one of the largest sources 
of uncertainty in our analysis, but it is also difficult to characterize. We also estimate an 
additional source of error in our selection of the emission volume that was chosen to represent 
one year worth of jet motion. Moreover, the correlation of in tegrated Br7 line lumin osity 
with stellar mass accretion luminosity has an intrinsic scatter (IMuzerolle et al.lll998bl ). and 
deviation in M out /M acc by a factor of several for a given target can result from using this 
method to derive the stellar mass accretion rate. Multiple observations of accretion indicators 
from large sample of T TS demonstrate the highly variable n ature of accretion activity in 
these young stars (e.g., IBary et al.l 120081 ; iNguyen et al.l 120091 ) . However, the effects from 
intrinsic time variation in the accretion and outflow properties are not an issue in our study 
because the mass outflow and mass accretion rates are derived simultaneously from the same 
data set. 

The blue-shifted jet emerging from DG Tau is arguably one of the best studied outflows 
associated with a youn g star. It was among the first cTTSs for which a collimated jet-like 



outflow was discovered (IMundt fc Friedlll983l ). Since its discovery, the HH 158 outflow from 



HST, ground-based adaptive optics systems, and spectro-imaging techniques 


L avalle v et al. 


1997; 


L avalle v- Fouq ue 


t et al. 


2000 


;] 


3acciotti et al. 


2000 


: Doueados et al.ll2002 


iBacciotti et al. 


2003; 


Pvo et al. 


2003; 


Takami et al. 


2004 


; 1 Coffey et al. 


2008). This outflow has observation- 



ally revealed the structure of collimated YSO jets in great spatial detail; they are typically 
comprised of an on-axis high velocity component (HVC) at radial velocities greater than 
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~50 km s 1 that can extend to spatial distances of hundreds of AU from the star, and an 
encompassin g shell of lower velocity gas that extends only about ~1 00 AU away from the 
central star faacciotti et al.lboool ; Ipvo et alJliooa iTakami et al.ll2004h . Although He 1 10830 
emission i s commonly a t tribu ted only to stellar winds in the inner ~10 AU environments 
of CTTS, ITakami et all ( 120021 ) detect He I 10830 emission from DG Tau at the jet velocity, 
spatially extended over 0.5 arcsec (~70AU) of the jet, thou gh their spectroastrom etry 
did not find exte nsion in the low- velocity He I 10830 emission. iBacciotti et all (120001 ) and 
Pyo et all (120031 ) described the blue-shifted HH 158 jet as having an "onion-skin" struc- 
ture, with high velocity low density gas on-axis, surrounded by successive layers of lower 
velocity, higher density gas. Integrated over the width of the jet axis, t he \ !,,„, I \ I,,, ,- values 



derived for the DG Tau jet from past stu dies lie in the range of 0.05-0.1 (jLavalley et al.lll997 



Bacciotti et all 120001 ; iCoffey et all 120081). though mos t of these estimat e s use large mass ac- 
cretion rate of ~10 6 M yr _1 ( IHartigan et al.lll995f ). IBacciotti et all ( 120031 ) analyzed the 



kinematics of DG Tau's jet and found evidence for velocity asymmetries across the jet axis, 
and they interpret this result as evidence for rotation within the inner jet channel. This 
controversial finding challenges the theories of magnetospheric interaction at the inner star 
+ disk + outflow boundary region. 

The outflow from DG Tau has served as a laboratory to test the physics of YSO jets, 
so it is fitting that our first discovery of spatially extended Br7 from cTTSs includes this 
system. In DG Tau, the spatially extended Br7 emission that we detect arises from the high 
velocity, on-axis gas in the blue-shifted outflow. Interestingly, this H I coincides very closely 
in projected l ocation with the re cently discovered X-ray emission in the jet (iGiidel et al. 
20051 . l2008h . iGiidel et all fj2008h found that the blue-shifted outflow from DG Tau has 
emission from soft X-rays out to spatially extended distances of several hundred AU from 
the central star. The temperatures of the soft X-ray emission regions associated with the 
extended X-ray jet are on the order of a few times 10 6 K. These temperatures pose a problem 
for understanding jet heating mechanisms. The standard jet shock models for YSO outflows 
cannot explain temperatures on the order of 10 6 K, and this argues strongly in favor of some 
form of non-thermal excitation of the gas; for example by Ohmic heating from magnetic 
dissipative currents (IGiidel et all 120081 ). The detected Br7 emission from the DG Tau blue- 
shifted outflow seems to correspond to gas nearby to these super-heated regions, and thus 
if strong magnetic fields significantly effect the physical conditions in the inner flows then 
our simple method of deriving M out assuming thermal Case B recombination of atomic 
hydrogen at temperatures of ~10 4 K would be correspondingly inaccurate. Further tests and 
comparison of the physical conditions in Br7 and X-ray emitting jets could help clarify the 
physics in the inner regions of the outflows. Unfortunately, as previously mentioned, HL Tau, 
HV Tau C and Haro 6-10 S are Class I protostars or cTTSs with circumstellar disks viewed 
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in a nearly edge-on configuration, and all three of these systems have very strong optical 
obscuration (see Table 2). It is likely not feasible to detect the soft X-ray emission from 
these jets to further test the physical conditions of the Br7 emission environments. 

6. Summary 

We have presented the results of our K-band integral field spectroscopy study of the 
Br7 line emission from eight T Tauri star systems. The key points of our study are: 

1) Using adaptive optics fed integral field spectroscopy, we have spatially resolved Br7 
line emission in the circumstellar environments around four of our eight survey target stars: 
DG Tau, Haro 6-10 S, HL Tau and HV Tau C. 

2) The spatially extended Br7 emission arises predominantly from the hydrogen re- 
combination regions associated with the inner Herbig-Haro outflows from these young stars. 
Only HL Tau shows a significant contribution of stellar Br7 emission scattered into our line 
of sight off of the outflow cavity walls. This emission from HL Tau has a similar morphology 
to the scattered light continuum flux, and it likely originated from the inner magnetospheric 
accretion region around the star. 

3) At some blue-shifted velocities, the spatially extended Br7 emission comprises ~20% 
of the detected Br7 (e.g., Haro 6-10). Although we spatially resolve Br7 emission from 
outflows in our high-contrast measurements, the majority of the integrated Br7 from most 
systems is spatially unresolved and may arise from the magnetospheric accretion processes 
at the location of the central stellar source (with the exception of HV Tau C). 

4) All of the Br7 emission that we detect above the continuum flux from HV Tau C is 
from the spatially extended emission, consistent with the location of the known Herbig-Haro 
outflow. HV Tau C is seen in continuum light as a scattered light edge-on disk nebulosity. 
The inner magnetospheric component of Br7 flux may be shielded from our line of sight by 
material in the inner edge of the circumstellar disk. 

5) Derivation of the stellar mass accretion rates from the relationship between Br7 line 
luminosity and mass accretion reveal M acc values that are typical of cTTSs. 

6) Detection of the spatially extended Br7 emission from the outflows in the DG Tau 
and Haro 6-10 S systems have allowed us to derive a value for the emitting hydrogen mass 
outflow rate using simple arguments applicable to hydrogen recombination regions. The 
corresponding values for M out /M acc that we derive are on the order of ~10-15%, consistent 
with many prediction from accretion- driven stellar winds and disk winds, while on the high 
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side in comparison to observationally derived mass outflow to accretion rate ratios from past 
studies. 

7) We find that in some young protostars, Br7 emission extended on spatial scales of 
greater than 0."1 (14 AU) can contribute ~ 10% of the flux to the detected integrated line 
emission (or more, as in the case of HL Tau and HV Tau C). 

Several NIFS datasets presented in this study were acquired during the early stages 
of instrument integration at Gemini North Observatory, and we are extremely grateful for 
the support of the NIFS teams at the Australian National University, Auspace, and Gemini 
Observatory for their tireless efforts during the instrument commissioning and system verifi- 
cation. This study is based on data from the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., on behalf of the international 
Gemini partnership of Argentina, Australia, Brazil, Canada, Chile, the United Kingdom, 
and the United States of America. 
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Table 1. Observing Log 



Star 


HH # 


Obs. Date 


PA of 


Exp Time/Coadds 


# Exp 


Total Exposure 


Note on 


Name 




(UT) 


Obs. 






Time 


Obs. 


T Tau 


HH 255 


2005 Oct 25 


0° 


5.3s/24 


36 


4580s 


AO Flexure Test 


RW Aur 


HH 229 


2005 Oct 22 


221° 


40s/l 


11 


440s 


AO Guide Test 


XZ Tau 


HH 152 


2005 Oct 25 


0° 


30s/l 


28 


820s 


AO+OIWFS* Guide Test 


DG Tau 


HH 158 


2005 Oct 26 


0° 


20s/6 


101 


12120s 


AO+OIWFS* Flexure Test 


HV Tau C 


HH 233 


2005 Oct 22 


114° 


900s/l 


3 


2700s 


System Sensitivity Test 


HL Tau 


HH 150 


2006 Feb 12 


146° 


900s/l 


3 


2700s 


0."2 Occulting Disk SV 


CW Tau 


HH 150 


2006 Dec 5 


146° 


900s/l 


3 


2700s 


NIFS PI GT 


Haro 6-10 


HH 150 


2007 Feb 7 & 8 


0° 


300s/l 


9 


2700s 


NIFS + LGS AO Queue 



"The NIFS On-Instrumcnt WaveFront Sensor (OIWFS) is used to correct spatial flexure in the observations. 
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Table 2. Stellar Parameters and Mass Accretion Rates 



Star 


Mass 


Luminosity 




A„ 


Integrated Br7 Flux 


Mace 


Mout 


Reference 




M© 


L 


K 


mag 


W/m 2 


M yr" 1 


M yr" 1 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


CW Tau 


1.1 


0.7 


4700 


2.2 


1.5xl0 -16 


1.3xl0 -08 






DG Tau 


2.2 


7.7 


4775 


3.3 


4.5xl0~ 16 


9.6X10" 08 


>1.2xl0-° 8 




Haro 6-10 S 


0.7 


1.8 


4000 


12.1 


7.2xl0~ 17 


6.7xl0~ 08 


>1.4xl0 -08 


Doppmann et al. 


HL Tau 


1.2 


3.0 


4400 


7.4 


2.1xl0- 16 


8.7xl0-° 8 






RW Aur A 


2.8 


12.9 


5000 


1.6 


7.2X10" 16 


1.4X10 -07 






RW Aur B 


1.2 


3.0 


4200 


1.6 


<5.1xl0" 18 


<4.1xl0- 10 






T Tau N 


2.1 


7.2 


5250 


1.5 


8.4xl0- ie 


1.5xl0-° 7 






XZ Tau A 


0.45 


0.4 


3400 


1.4 


<7.7xl0" 18 


<1.9xl0 -09 






XZ Tau B 


















Haro 6-10 N 


0.6 


0.8 


4100 


29.0 


1.7xl0 -17 


4.6xl0" 08 




Doppmann et al. 
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Fig. 1. — (a) The 2.16 /im continuum flux level from DG Tau with contours of [Fe II] 
emission overplotted, designating the outflow location, (b) the continuum subtracted point- 
source Br7 flux and (c) the continuum and point-source subtracted map of spatially extended 
Br7 from DG Tau. Also included here is a continuum and point-source subtracted map of 
spatially extended Br7, integrated over only three blue-shifted velocity channels (d). The (d) 
panels show that much of the spatially extended Br7 emission is blue-shifted. The locations 
designated as "B" and "C" had 1-D spectral traces extracted, these are presented in Figure 
5. 
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Fig. 2. — (a) The 2.16 /im continuum flux level from Haro 6-10 with contours of [Fe II] 
emission overplotted, designating the outflow location, (b) the continuum subtracted point- 
source Br7 flux and (c) the continuum and point-source subtracted map of spatially extended 
Br7 from Haro 6-10 S. Also included here is a continuum and point-source subtracted map of 
spatially extended Br7, integrated over only three blue-shifted velocity channels (d). The (d) 
panels show that much of the spatially extended Br7 emission is blue-shifted. The locations 
designated as "B" and "C" had 1-D spectral traces extracted, these are presented in Figure 
5. 
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Fig. 3. — (a) The 2.16 /im continuum flux level from HL Tau with the 0."2 occulting disk 
in the beam with contours of [Fe II] emission overplotted, designating the outflow location, 
(b) the continuum subtracted point-source Br7 flux and (c) the continuum and (d) point- 
source subtracted map of spatially extended Br7 from HL Tau. Also included here is a 
continuum and point-source subtracted map of spatially extended Br7, integrated over only 
two blue-shifted velocity channels (d). The (d) panel shows weak blue-shifted Br7 emission 
that arises from the same spatial location as the [Fe II] emission from the outflow. The 
locations designated as "B" and "C" had 1-D spectral traces extracted, these are presented 
in Figure 5. 
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Fig. 4. — (a) The 2.16 fim continuum flux level from HV Tau C with contours of [Fe II] 
emission overplotted, designating the outflow location, and (b) the continuum subtracted 
point-source Br7 flux. The locations designated as "B" and "C" had 1-D spectral traces 
extracted, these are presented in Figure 5. 
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Fig. 5. — Br7 flux versus velocity plots at the continuum peak locations and "Position B" 
and "Position C" (from plots 1-4) for DG Tau, Haro 6-10 S, HL Tau and HV Tau C. For 
DG Tau, HL Tau, and HV Tau C, the more spatially extended Br7 emission measured at 
the "Position C" location is fainter and more blue-shifted than the less extended flux found 
at "Position B". For Haro 6-10 S, the "Position C" flux measures a bright knot in the 
Herbig-Haro flow and is stronger and more blue-shifted than the "Position B" emission. 
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Fig. 6. — The standard deviation of the continuum subtracted Br7 flux plotted versus 
distance from the central continuum peak for HL Tau (a), DG Tau (b), Haro 6-10 S (c) 
and HV Tau C (d). The overplotted dashed line represents the level of spatially extended 
continuum and point-source subtracted Br7 flux. The ratio of the spatially extended flux 
to the standard deviation of point-source subtracted emission (dashed line to solid line) 
represents a measure of the signal-to-noise of the spatially extended Br7. The peak of the 
S/N for extended Br7 is ~11 for HL Tau, ~6 for DG Tau, ~25 for Haro 6-10 S, and ~11 
for HV Tau C. 
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Fig. 7. — The standard deviation of the continuum subtracted Hi" (lux plotted versus 
distance from the central continuum peak for CW Tau (a), T Tau (b), XZ Tau (c) and RW 
Aur (d). The overplotted dashed line shows the signal from point-source subtracted Br7 
flux with increasing distance from the central star (e.g., XZ Tau A, RW Aur A, and T Tau 
N, in the case of the multiples). Flux peaks at 0."65 and 0."25 seen in T Tau and XZ Tau 
represent the positions of the stellar companions in these systems. The ratio of the dashed 
line to the solid line provides a measure of signal-to-noise in spatially extended Br7 flux, 
though the companions complicate the detection. It is clear from these plots that XZ Tau B 
and T Tau S have significant Br7 emission above their continuum flux. In these four cases, 
the S/N of spatially extended Br7 is less than ~3 and there is no convincing evidence for 
spatially extended emission in any of these TTSs. 



